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Abstract 

Recently, an exact Green's function of the diffusion equation for a pair of 
spherical interacting particles in two dimensions subject to a backreaction bound- 
ary condition was derived. Here, we use the obtained Green's function to calcu- 
late exact expressions for the survival probability, the time-dependent reaction 
rate coefficient for the initially unbound pair and the survival probability of the 
bound state in the time domain. Moreover, we derive an exact expression for 
the off-rate. 

1 Introduction 

In [T], an exact Green's function (GF) g{r,t\ro) of the two dimensional (2D) 
diffusion equation for a pair of spherical interacting particles was derived. The 
GF satisfies the backreaction boundary condition [2l [3l |4j [1] 

d 

2TraD—g{r,t\ro)\r=a = Kag{a,t\ro) - - S{t\ro)]. (1.1) 
or 

Here, D is the sum of the particles' diffusion constants, r is the inter-particle 
distance and a is the encounter radius. Furthermore, Ha and Hd denote the 
intrinsic association and dissociation rate constants, respectively. Finally, S{t\ro) 
refers to the probability that a pair of molecules with initial distance vq survives 
until time t. 



Sit\ro) = 27r / g{r,t\ro)rdr (1.2) 



= 1 - 



27raZ?^ ^g{r,t'\rn)\r=adt' . (1.3) 
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It was shown that the associated GF takes the form [I] 

9ir, t\ro) = — / e-^"'*T(x, r)T{x, tq) x dx, (1.4) 



27r .„ 

where we introduced the funtions 

rj., ^ _ Jo{rx)l3{x) ~YQ{rx)a{x) 

[a(a;)2+/3(a;)2]i/2 ' ^''^^ 

and 

a{x) :— {x'^ — K]j)Ji{xa) + hxJo{xa), (1-6) 
I3{x) := {x'^ ~ KD)Yi{xa) + hxYoixa). (1.7) 

Jq, Ji,Yo,Yi denote the Bessel functions of first and second kind and of zeroth 
and first order, respectively Furthermore, by definition we have 

KB := (1,9) 

Knowing the GF allows us to derive further important quantities, notably, the 
survival probability. For the following calculations, it turns out to be convenient 
to introduce the function 

Pfa r) ^lAr(x r) ^ M^mx)-Y,{xr)a{x) 

xdr [a(a;)2+/3(a;)2]i/2 ' ^'■''^> 

where we used [5] 

Joix) -^Mx) = -Mx) (1.11) 

Y^{x) A.Yo{x) = -Y,{x). (1.12) 

The survival probability may be calculated according to (jl.2p or, alternatively, 
using (jl.3|) . Either way, we find 

/>oo 

S{t\ro) = l-a e-^*^'P(a;,a)T(x,ro)da;. (1.13) 
Jo 

To obtain (|1.13D via (jl.3p we used the integral identity 

P{x,a)T{x,ro)dx = 0. (1.14) 



p.l4p can be derived by substituting the GF (|1.4p into the boundary condition 
(fTTj) . Then, using ([ITU)) . (fTTT]) and ([TT^ . as well as [S] 

Y^{x)Jo{x) - 4ix)Yo{x) = — , (1.15) 

TTX 

we are directly led to (|1.14p . From (jl.l3l) we can easily conclude that the re- 
versible survival probability approaches unity for large times 

lim S{t\ro) = 1, (1.16) 
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implying that in 2D the uhimate fate of an isolated pair for the reversible reaction 
is always dissociation, as in the ID and 3D case [2l[3j|4]. 

Turning now our attention to the initially bound pair, we use the notation 
* to indicate the bound state and we let g{r,t\*) denote the Green's function 
describing an initially bound pair that is found separated by a distance r at a 
later time t, cp [21I2- Knowledge of this GF is desirable, because it enables us 
to compute the average lifetime of the bound state and, hence, the off-rate, as 
we shall see later. 

We note that g{r,t\*) is related to the GF g{*,t\rQ), which describes an 
initially unbound pair that is bound at a later time t, by the detailed balance 
condition [31 13] 

i^dg{*,t\r) = Kag{r,t\*). (1-17) 

Furthermore, we have 

g{*,t\ro) + S{t\ro)^l. (1.18) 

(jl.lSp makes it evident that, technically speaking, the GF g{*,t\ro) represents a 
probability rather than a probability density. 
From ([113)) . ([LTf)) and ([LTS)) it follows that 

g(r,t\*)^—a e^°'^'P{x,a)T{x,r)dx. (1.19) 

'^a Jo 

Upon direct integration we arrive at the probability 5(^1*) that an initially bound 
pair is unbound at time t > [3 H] 

S{t\*) = 2TT rg{r,t\*)dr = l-2TT—a^ e-^'''' P\x, a)-dx. (1.20) 

Ja JO 

Note that S'(0|*) = 0, due to ([X^ . 

In [2 it was demonstrated that S{t\*) relates to S{t\rQ) in the time domain 
in the following way: 

S{t\*)^Kdf e-'""*' S{t-t'\a)dt'. (1.21) 
Jo 

Using ([T3)) . ([TTTU)) . ([TTT3)) . ([TTT5)) . and ([L^ . we can now exphcitly show that 
(|1.2ip is indeed satisfied in 2D. 



2 Rates 

In the Smoluchowski approach, steady-state rate constants can be calculated in 
two related ways. One method employs the irreversible time-dependent solution 
of the diffusion equation to calculate the time-dependent rate coefficient ki„(t). 
The steady-state rate constant fc?" is then given by the long time limit of kirr{t) 

k°}l := lim k,,{t) (2.1) 

t—^oo 

While in 3D this procedure yields the classic result fc™ = inaD (for purely 
absorbing boundary conditions), in 2D one obtains the result that kf^ actually 
vanishes. This issue is also reflected in the second approach, which is based on 
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steady-state solutions of the diffusion equation. In 2D, however, the steady-state 
solutions are logarithmic, and hence, no boundary condition can be imposed at 
infinity. As a remedy, several procedures have been proposed O [71 [HI IS] ■ All of 
them give rise to another length scale 6 > a in addition to the encounter radius. 
Furthermore, they lead, up to a constant, to the same expression for the on-rate 
constant 

1 1 \xi{h/a) + 5 



Using iCcq :— — K^^k^^ one arrives at an expression for the off-rate 



(2.2) 



fcoff Kd '27tD 



Here 5 denotes the above mentioned constant that assumes a varying value de- 
pending on the chosen method. 

In the following we will use the reversible GF (|1.4p to address these issues. 
The time-dependent reaction rate coefficient can be obtained from the survival 
probability S{t\ro) [3 [J 

k{t) = 2TTaD-^^S{t\ro)\ro=a- (2.4) 
Using f|1.13p . we obtain the exact expression in the time domain 



oo 



k{t)^2Tra^D e-^^'-' P\x,a)xdx. (2.5) 



It follows immediately from (|2.5I) that the long time limit of k{t) vanishes. 

As shown in the appendix, the integral P^{x,a)^ can be computed an- 
alytically and we find 

P2(.,a)^ = ^^^. (2.6) 

With the help of this equation and (|2.5p we recover the correct expression for 
the equilibrium constant K^q 

/ kit)dt = ^ =: /^eq. (2.7) 

Alternatively, identity (|2.7p may be deduced in the following way. According 
to reference T the time-dependent reaction rate can also be calculated by using 
the survival probability ^(tl*), cp (|2.4p 

dS{t\*) 

at 



m = K^q^^. (2.8) 



Now, given (ITTTO)) and (IL^ . one can explicitly show that (|^ and (|^ yield 
the same result. Then, taking into account that S'(oo|*) = 1 and ^(Ol*) = 0, 
(|^ is immediately implied by (|^ . 

The steady-state dissociation rate constant, or off-rate, is defined by [3] [^ 

fcoff -r-i, (2.9) 
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where r refers to the average hfetime of the bound state, r can be calculated 
according to [3j |4] 



/■oo 

/ [I - S{t\*)]dt. (2.10) 
Jo 



Using (|1.20p and after transforming from the dummy integration variable x 
to the dimensionless integration variable ^ :— ax we arrive at 

..2,-^1 fM*. ,2.11, 



D Jo £_ 

where we have introduced the function 

f(0 ■■= (2.12) 

for notational convenience. We note that the integrand in (|2.1ip depends on 
the dimensionless constants h = ha, ku :— k^o^ only, due to the transformation 
a; — > ^ = xa. Since 

lim/(O = ^^^0, (2.13) 

it follows that the integrand is singular at the lower endpoint ^ = of the 
integration interval, which implies that the integral does not exist. However, one 
can still assign a well-defined value to integrals of that type upon regularizing 
them in the sense of Hadamard's one-sided finite-part integrals [10] 



By definition, the following finite part integral yields ^Oj 



(2.14) 



Inc (2.15) 

We would like to point out that one has to split the integral at oo > ^ = c > in 
(j2.14p , because otherwise the regularization would introduce another singularity 
due to the occurrence of In oo terms. It is instructive to compare the resulting 
expression for r with (j2.3p . Clearly, the third term on the rhs of (j2.14p gives rise 
to precisely the logarithmic contribution if one set c = h/a. Moreover, one can 
show numerically that the other two integrals in (|2.14p yield 

-^d^+ d^^^ + C{c)-r^, (2.16) 

where C(c) depends on the choice of c only. Taken together this means that 
one seemingly recovers the result (|2.3p . However, it is easily demonstrated that 
the regularized integral is actually independent of the choice of c. Thus, in 
particular, one can choose c = 1, which shows that the logarithmic contribution 
actually vanishes and one obtains 



— - — -f if — 



Ke^T^ZT^, (2-17) 



where C := C{1) « 0.11593 We find the result that the off-rate does not 

depend on the encounter radius. 
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3 Appendix 



In this appendix we will compute the integrals that are needed in the main text. 
The GF has to satisfy the initial condition [I] 

5(M = 0M = ^^, (3.1) 
i.e. in the case considered here we have 

l^fTi.,r)Ti.,ro).d.= '-^. (3.2) 

From p.2p it follows by direct integration over dr'r' 

' P{x,r)T{x,ro)dx^(^lu^ l^^'^ (3.3) 

where we have used (jl.l4p . Furthermore, to perform the integral we used [5] 

-^[a;''J,(x)] (3.4) 
ax 

and 

■^[x^Y^ix)] ^ x^Y.^i. (3.5) 
ax 

Next, we apply the derivative to the boundary condition (jl.ip and subse- 
quently integrate from to t. Taking the limit t ^ oo leads to 

P(x,a)P(x,ro)xdx — — / T(x,a)P(x,ro)dx + 

2TTaD Jo 

^ r P{x,a)P{x,ro)-. (3.6) 
D Jo X 

The integral on the Ihs has to vanish, as can be seen by taking into account 
(fTTl) and 

/ P{x,a)P{x,ro)xdx — — / P{x,a)T{x,ro)dx. (3.7) 
Jo clfo Jo 

Furthermore, we just calculated the first integral on the rhs, cp Eq. p. 31) . Hence, 
we arrive at ^ 

P(.,a)P(.,.o)- = -^. (3.8) 
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